Thermodynamic equilibrium conditions of graphene films on SiC 
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First-principles surface phase diagrams reveal that epitaxial monolayer graphene films on the Si 
side of 3C-SiC(lll) can exist as thermodynamically stable phases in a narrow range of experimentally 
controllable conditions, defining a path to the highest- quality graphene films. Our calculations 
are based on a van der Waals corrected density functional. The full, experimentally observed 
(6v^ X 6a/3)-R30° supercells for zero- to trilayer graphene are essential to describe the correct 
interface geometries and the relative stability of surface phases and possible defects. 



The growth of wafer-size graphene films on a semicon- 
ducting substrate is a first step towards graphene based 
electronics. The semiconductor SiC as a substrate may 
hold the key to device applications. Here, well-ordered 
graphene films can be grown directly on a semiconduct- 
ing substrate by a simple process (Si sublimation from 
the surface, e.g., Refs. [iMzl) and the standard tools of 
semiconductor technology can be used for further ma- 
nipulation. Indeed, graphene-based devices and even in- 
tegrated circuits [8|, |9| were already created on the Si 
side of SiC substrates. However, controlling the precise 
thickness of graphene films is important to minimize the 
coexistence of monolayer graphene (MLG) and bilayer 
graphene (BLG) . [SL IsL lloLll l| MLG areas exhibit no band 
gap, while BLG areas do.j^ A recent, joint experimental- 
theoretical study finds particularly high local resistances 
across monolayer-bilayer graphene junctions on the same 
surface, [12| a possible contributing factor to low carrier 
mobilities in graphene on SiC(0001).[i,[i2| 

The growth of graphene on SiC is special in the sense 
that, instead of offering one or more of the components 
from the gas phase, graphene areas are formed by con- 
trolled sublimation of Si from the surface. [ll,l2| Graphene 
films grown under ultrahigh vacuum conditions are typ- 
ically inhomogeneous.jsl, l5[ To improve their quality is 
therefore a major and ongoing experimental goal.jsl |5| 
In the past, th e appea rance of different phases was often 
interpreted [5|, [isl, [l^ as intuitive, successive intermedi- 
ates formed by an outgoing Si flux that ultimately leads 
to bulk-like graphite layers. If graphene films of various 
thicknesses were a result of a purely kinetically limited 
Si sublimation process (controlled by growth tempera- 
ture and time), improving upon the homogeneity of the 
layer thickness would be difficult. It would be helpful 
if there were a set of thermodynamic equilibrium condi- 
tions (e.g., temperature T and partial pressures p of Si 
and C [l5|) at which a desired film thickness were ther- 
modynamically stable over all others. As long as T is 
sufficiently high to overcome the relevant kinetic barri- 
ers, large-scale ordered films of that particular thickness 
could then be formed simply by finding and maintaining 
(p, T) near these equilibrium conditions. 

In fact, experimental evidence is mounting that con- 
trolled (p, T) conditions do indeed aid the growth of 



graphene on SiC. Ref. [16| demonstrates that the carbon- 
rich "zero-layer graphene" (ZLG) or "buffer layer" pre- 
cursor phase [1, 2, [isl, [l7| (not yet graphene) on the 
Si face is a reversible thermodynamic equilibrium phase 
at high T with a controlled disilane background pres- 
sure. Reversibility is much harder to demonstrate once a 
complete graphene plane is formed, [18| but an increased 
growth temperature in an Ar background buffer gas does 
lead to much improved MLG film homogeneity. [3| Excel- 
lent wafer-size MLG films are also reported for growth in 
a confined cavity that may retain a finite, well-defined Si 
background pressure as Si evaporates from the surface. [5| 
Finally, a well-defined graphene precursor phase on the 
C face at finite disilane background pressure was re- 
ported very recently. [6] What is still not clear, however, is 
whether MLG itself is an equilibrium phase under certain 
conditions. If so, one could ideally facilitate the growth 
of MLG but not BLG on SiC(lll). 

We here present first-principles theoretical evidence 
that such equilibrium conditions indeed exist for MLG 
(and, possibly, even BLG) on the Si face of SiC (111). 
We employ density-functional theory (DFT) using the 
van der Waals (vdW) corrected jl9| PBE density 
functional [20] (called PBE+vdW throughout this work). 
Strategies to include vdW effects have become an ac- 
tive scientific area of their own (see, e.g., Refs. |2lM26J 
and many references therein). However, for the large, 
complex carbon-rich interfaces of interest here, no un- 
ambiguously quantitatively improved approach over the 
level of theory used here exists to our knowledge, i.e., 
the present incorporation of vdW terms reflects the state 
of the art. For the structure of the relevant bulk phases, 
the impact of different standard functionals is well under- 
stood and systematic. Predicted lattice parameters (see 
supporting information (SI) for full reference data|27J) 
for 3C SiC are within 1 % of PBE+vdW and exper- 
iment: aac-sic =4.33 A, 4.38 A, 4.36 A, respectively, 
for the local density approximation (LDA), PBE with- 
out vdW, and PBE+vdW. Similary small discrepancies 
arise for diamond C and for the in-plane lattice parame- 
ter of graphite: agraphite=2.45 A (LDA), 2.47 A (PBE), 
2.46 A (PBE+vdW). Zero-point corrections (ZPC) to the 
lattice parameters are also below 1 %.[27|. For the in- 
terplanar lattice parameter c of graphite, vdW effects 



must be included into the PBE functional (LDA: 6.65 A, 
PBE+vdW: 6.66 A, but 8.65 A for PBE). On a tech- 
nical level, our calculations are based on the FHI-aims 
all-electron code[28|, |29| with "tight" numerical settings 
and the massively parallel ELPA eigensolver library [30| 
to guarantee accurate total energies for the very large 
structure sizes involved (see SI for details [27]). 

In experiment, one encounters a series of phases on the 
Si side of 3C-SiC(lll) (experimentally also observed on 
the Si side of 4H-SiC(0001) and 6H-SiC(0001)) when go- 
ing from a Si-rich to a C-rich regime. In UHV, a (3x3) 
Si-rich layer [31] can be prepared. Upon Si sublima- 
tion, a simpler (>/3 x \/3)-R30° bulk-terminated surface 
with one adsorbed Si adatom per three unit cells fol- 
lows. [32I, |33| Removing yet more Si creates a partially 
covalently bonded carbon "buffer layer," the ZLG phase, 
with a large, commensurate unit cell: One (13x13) 
honeycomb graphene-like supercell (338 C atoms) on a 
(6V^x6V3)-R30° mesh (108 Si and 108 C atoms per bi- 
layer) of the underlying SiC substrate, jl^ Compared to 
a graphene plane in graphite, the lattice match is al- 
most strain-free (experiment: 0.2 % at T=0 K [3J, |35|, 
PBE+vdW: 0.1 % [13). The ZLG phase does not 
yet exhibit the electronic properties of actual freestand- 
ing graphene. [13] Further heating detaches the ZLG C 
plane from the substrate to form MLG and a new C- 
rich layer underneath. [36|, |37| The same process can be 
continued to successively form BLG [38| and multilayer 
graphene films. Importantly, our calculations address 
the graphene-like films in their experimentally observed, 
large commensurate (6\/3 x 6>/3)-R30° supercells, us- 
ing slabs containing six Si-C bilayers under each recon- 
structed phase (1,742 up to 2,756 atoms for ZLG up to 
three-layer graphene (3LG), respectively). The top three 
SiC bilayers and all planes above are fully relaxed (resid- 
ual energy gradients: 8-10~^ eV/A or below). 

Figures [T^-c show the ZLG, MLG and BLG phases to- 
gether with key geometry parameters predicted at the 
level of PBE+vdW. Since all planes are corrugated. |39l- 
4l| histograms for the atomic z coordinates are included. 
The interface geometry stays essentially the same despite 
the addition of more graphene planes. In the MLG phase, 
we see a significant buckling of the topmost graphene 
layer {^OAl A between top and bottom of the plane). 
This strong buckli ng is qu alitatively consistent with ex- 
isting STM images (iJI-li^. In the BLG phase, the corru- 
gation is slightly reduced, but the two top planes are still 
buckled by 0.24 A and 0.32 A, respectively. This buckling 
reflects some coupling to the covalently bonded interface 
C-rich plane, which is much more corrugated (^0.8 A 
in our work, similar to experimental estimates 451 146|). 
The observed graphene interplanar distances near the in- 
terface are slightly expanded compared to experimental 
bulk graphite (3.34 A [35]) and in good qualitative agree- 
ment with estimates from scanning tunneling microscopy 
|44| and transmission electron microscopy. |47|] 
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FIG. 1: a-c: Geometry and key geometric parameters for 
the three phases ZLG, MLG and BLG on the Si face of 3C- 
SiC(lll), as determined by DFT-PBE+vdW. The histograms 
distribute the atoms at each z coordinate (relative to the 
fourth bilayer of the substrate) into bins of 0.01 A width, 
giving an impression of the overall corrugation in each plane. 



We have also compared our findings to geometries 
for the straight PBE functional (no vdW correction), 
and for the local-density approximation (LDA). In PBE, 
the C-G interplanar distances are unphysically expanded 
(4.42 A for MLG). In contrast, the LDA geometry of 
the carbon planes agrees qualitatively with PBE+vdW, 
although LDA incorporates no long-range vdW interac- 
tions. The first qualitative geometry difference between 
the PBE+vdW and LDA treatments appears in the Si 
part of the top Si bilayer, where some Si atoms (those 
with dangling bonds pushing against the 7r-bonded parts 
of the C interface plane [48]) are located much deeper (by 
?^0.3 A) in PBE and PBE+vdW than in LDA. A direct 
crystallographic verification would be desirable. 

In a grand canonical formalism, the possible equilib- 
rium conditions for different surface phases can be rep- 



resented by the chemical potentials of C and Si, /ic and 
/isi- In experiment, /ic and /isi can be manipulated, e.g., 
through the substrate temperature and background pres- 
sure of gases that supply Si or C.[5|, |6|, [iSl, ll6|. Consider 
the surface energies 7 of a two-dimensional periodic SiC 
slab with a C face and a Si face. In the limit of sufficiently 
thick slabs, we have: 



7Si-face + 7C-face = -J (^' 



slab 



A^Si/iSi - A^c/ic) . (1) 



Nsi and A^c denote the number of Si and C atoms in the 
slab, respectively, and A is the chosen area. All our sur- 
face energies are in eV per area of a (1x1) SiC unit cell. 
The letter E denotes total energies for a given atomic ge- 
ometry throughout this work. In our actual calculations, 
we always choose a fixed H-terminated C-face geometry, 
which cancels out for all surface energy differences related 
to the Si face. 

The major experimental (T,p) dependence during 
growth arises through the reservoirs of Si and C, which 
define /isi and /ic • [ISl, ll6| Thus, a precise control of back- 
ground gases as reservoirs (for instance, disilane |16| ) 
is desirable, even if calibration variations [49| may re- 
quire exact (T,p) ranges to be adjusted separately for 
a given growth chamber. The actual growth process pro- 
ceeds by Si out-diffusion from underneath already formed 
graphene planes. Yet, the external Si reservoir back- 
ground pressure still matters in equilibrium: As long as 
the diffusion path to the outside remains open, so does 
the inward diffusion path, and near equilibrium with the 
reservoir gas can be achieved. During intermediate stages 
of the formation of a new graphene plane, [47] such diffu- 
sion paths must be available. 

In principle, we could further include the much smaller 
(p, T) dependence of the solid phases by focusing on 
Gibbs free energies G{T^p) in the (quasi-) harmonic ap- 
proximation instead of E. However, quantifying this 
T dependence precisely would here necessitate accurate 
phonon calculations for structure sizes of the order of 
^2,000 atoms, a task that is computationally prohibitive 
at present. ZPC are small for the bulk phases (see SI 
[27|). Still, the possible small contribution of finite T 
stresses at the growth conditions [50| is kept in mind when 
interpreting our calculated results below. 

In equilibrium with a stable SiC bulk, /ic and /isi are 
linked through 
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TTibulk 



TTibulk 



2Ai7/(SiC). (2) 



The energies are per atom, and AH f {SiC) is the forma- 
tion enthalpy of SiC with respect to the elemental C and 
Si. The bulk phases define the chemical potential limits 
within which the SiC crystal is stable against decompo- 
sition into bulk Si or C: /ic < E^""^^ and /isi < E^^"^^, 
leading to 



ZTibulk 



AHfiSiC) < lie < E^ 



bulk 



(3) 



and analogous for Si. The diamond structure for Si is 
the appropriate bulk phase, but for C, there is a close 
competition between diamond and graphite. [SjJ, l52| We 
thus include both phases in our analysis. 

The experimentally reported energy difference between 
diamond and graphite at T=0 K is 25 meV/atom [5l|. 
Based on the potential energy minima (no ZPC) graphite 
is found to be more stable than diamond in PBE+vdW 
by 60 meV/atom (Fig. [2^). This is qualitatively con- 
sistent with the extrapolated experimental phase hi- 
erarchy. In plain PEE graphite is overstabilized by 
130 meV/atom. In LDA, both phases are similarly sta- 
ble: Considering only the potential energy surface, dia- 
mond is slightly more stable (by 12 meV), but already 
the inclusion of ZPC[52| would neutralize this balance 
(graphite more stable by 3 meV/atom[53|). 

The surface energies of the known surface phases of 
SiC (Si face) are shown as a function of A/ic = jic — 
£;buik -^ p-g [2^ £q^ PBE+vdW. The most stable phase 
for a given value of A/ic is that with the lowest surface 
energy. Going through Fig. [2^ from left to right, we 
find the expected broad ranges of stability for the Si-rich 
(3x3) and (V^ x V^)-R30° phases. Just before the C- 
rich limit (bulk graphite) is reached, there is a crossover 
first to ZLG, then to MLG, and even to a very narrow 
slice of the BLG phase. As an additional bound, ABA- 
stacked trilayer graphene (3LG) is also shown, crossing 
BLG within 1 meV of bulk graphite. 

While the respective stability ranges are narrow (in- 
set of Fig. [2^: 4 meV, 5 meV and <1 meV for ZLG, 
MLG and BLG, respectively at the chemical potential 
axis), but it is important to recall that narrow chemical 
potential ranges do not necessarily correspond to narrow 
experimental conditions: A drastic change in the number 
of Si (A/'si) and C (A^c) atoms can correspond to a small 
change of /i. For instance, one would first have to remove 
all Si from the SiC crystal to cross beyond the graphite 
stability line in equilibrium. However, the surface-energy 
differences between the different phases are also rather 
small (a few tens of meV per (1x1) SiC surface unit 
cell). The primary approximations that we cannot sys- 
tematically improve in our calculations are the density 
functional used, as well as possible small temperature- 
dependent surface strain effects (see above). The key 
message of Fig. [2^ is thus that MLG and its related 
phases all appear at least as very near equilibrium phases, 
a fact that is nonetheless critical for a qualitatively cor- 
rect understanding of their growth and properties. 

What we can do is to show how our results would 
be affected by different density- functional treatments. 
We have thus recomputed the surface phase diagram up 
to the MLG phase for two widely used functionals in 
Figs. [2)3 and c: the plain PBE functional, which lacks 
long-range vdW interactions and should thus yield un- 
trustworthy results, and for the often used LDA func- 
tional. As expected, the absence of vdW tails in the plain 
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FIG. 2: Comparison of the surface energies for five different reconstructions of the 3C-SiC(lll) Si side, relative to the bulk- 
terminated (1x1) phase (always unstable), as a function of the C chemical potential within the allowed ranges (given by 
diamond Si, diamond C or graphite C, respectively), (a) PBE+vdW, (b) PBE, (c) LDA. The shaded areas indicate chemical 
potential values outside the strict thermodynamic stability limits of Eq. [3] 



PBE functional changes the phase diagram drastically. 
Due to the overstabilization of graphite (130 meV/atom), 
its stability line moves significantly further to the left, as 
does the crossover point between ZLG and MLG. As a 
result, neither ZLG nor MLG becomes stable over the 
competing Si-rich (V^ x V^)-K30° phase in PBE, in 
outright contradiction to experiment. [16] In the LDA- 
derived phase diagram, the most significant change com- 
pared to PBE+vdW is the apparent incorrect stability 
hierarchy of graphite vs. diamond (without ZPC). If the 
diamond line were discounted, our calculations show the 
ZLG-MLG crossover point almost exactly on the graphite 
line. Still, even taking LDA at face value implies the ex- 
istence of T-p conditions very close to equilibrium for 
MLG, making the experimental search for such condi- 
tions promising. 

Figure [2] thus shows the most important point of our 
paper: The existence of equilibrium or near-equilibrium 
chemical potential ranges for ZLG, MLG, and even BLG, 
corresponding to specific T/p conditions in experiment. 
For each phase, this finding proves the potential for much 
better growth control than what could be expected if each 
phase were just a necessary (but not thermodynamically 
stable) kinetic intermediate. While true reversibility for 
actual MLG may be hard to achieve [16|, ll8| (the reverse 
growth process, disassembling a fully formed graphene 
plane would be kinetically difficult), the active forward 
growth process from MLG to BLG under Si out-diffusion 
should still be limitable by appropriate T/p conditions. 
A macroscopically homogeneous surface very close to 
pure-phase MLG should thus be achievable in principle. 

Figure [2] shows unambiguously the importance of a 
consistently accurate numerical treatment of the exper- 
imentally observed phases in their large unit cells. It 
would obviously be much more economical to consider 
smaller-cell approximant phases to the true (6\/3x6\/3)- 
R30° supercells. However, the residual artificial strain 



and inadequate bonding in those phases are too large 
for meaningful surface energy comparison s. |54l For in- 
stance, the popular (\/3 x >/3)-R30° [48|, l55| approxi- 
mant would intersect the graphite stability line at a sur- 
face energy 7! E;buik =0.15 eV, far above the actu- 

/^G graphite 

ally stable phases. Likewise, a slightly rotated (5x5) 
approximant to the ZLG phasel56| (a periodicity some- 
times seen in experiment! 171. |57|) would intersect at 
(7l/7^=£;buik =-0.35 eV), still higher by 0.06 eV than 

H'C graphite 

even the closest competing Si-rich phase, the (>/3 x VS)- 
R30° Si adatom phase (7L^^£;buik = -0.41 eV). The 

/J-C graphite 

(5x5) phase is either a nonequilibrium phase, or its struc- 
ture is not the same as that assumed in Ref. |56| . 

The true problem with artificially strained approxi- 
mant phases is that the resulting strain can obscure other 
electronically relevant properties, such as the energet- 
ics of defects. As an example, we consider a specific 
class of C-rich defects suggested as an equilibrium fea- 
ture of the ZLG phase in Ref. [58]. Two different de- 
fect positions, "hollow" and "top" were suggested. [58| 
Indeed, both would be more stable than the hypothetical 
{VS X \/3)-R30° ZLG approximant when included there 
in a (3 X 3) arrangement as done in Ref. [58|: —1.75 eV 
per defect for the hollow position, —2.93 eV per defect for 
top, both at /ic = ^gmphite- However, the same defects 
are unstable when included into and compared to the cor- 
rect (6V^ X 6\/3)-R30° ZLG phase: +5.28 eV per defect 
for hollow, +5.27 eV for top, again at /ic 
SIj27| for structure and other details). 



^gmphite (^^^ 



In conclusion, we can now rationalize some specific 
growth-related observations: 

(1) When simply heating a sample in UHV, the back- 
ground pressures of Si and C are low and ill-defined. The 
observed inferior morphologies and wide variations of ex- 

I — I r"^ 

perimental conditions [49|,|59| are consistent with this pic- 
ture. 



(2) Much more homogeneous growth can be achieved in 
an Ar atmosphere, [3] although MLG/BLG phase areas 
still coexist. A uniform background partial pressure of 
Si, however, will not be strictly guaranteed. 

(3) The observed thermodynamic ZLG stability and im- 
proved growth of MLG by controlling a disilane reservoir 
[16| is fully consistent with our findings. 

(4) The use of a confined cavity to control the Si flux away 
from the sample reportedly yields excellent wafer-size 
films. [5[ Maintaining a controlled Si partial pressure at 
constant temperature is most likely the important step. 

(5) Finally, similar considerations may also aid the much 
more difficult growth of MLG on the C face.[5|,[60| While 
our work is restricted to the Si face, a well-defined in- 
terface layer on the C face at finite disilane background 
pressure was proposed very recently. [6| This finding is an 
excellent additional indication that near-equilibrium sur- 
face conditions are indeed the key to the best possible 
epitaxial growth of graphene on SiC. 



[1] 

[2] 
[3] 

[4] 
[5] 

[6] 
[7] 
[8] 

[9 

[lo; 
[11 

[12 
[13 
[14 
[15 
[16 
[17 



A. Van Bommel, J. Crombeen, and A. Van Tooren, Surf. 
Sci. 48, 463 (1975). 

I. Forbeaux, J.-M. Themlin, and J.-M. Debever, Phys. 
Rev. B 58, 16396 (1998). 

K. V. Emtsev, A. Bostwick, K. Horn, J. Jobst, G. L. Kel- 
logg, L. Ley, J. L. McChesney, T. Ohta, S. A. Reshanov, 
J. Rohrl, et al., Nature Materials 8, 203 (2009). 
C. Riedl, C. Coletti, and U. Starke, J. Phys. D: Appl. 
Phys. 43, 374009 (2010). 

W. A. de Heer, C. Berger, M. Ruan, M. Sprinkle, X. Li, 
Y. Hu, B. Zhang, J. Hankinson, and E. Conrad, Proc. 
Nat. Acad. Sci. 108, 16900 (2011). 

N. Srivastava, G. He, Luxmi, and R. M. Feenstra, Phys. 
Rev. B 85, 041404 (2012). 

G. R. Yazdi, R. Vasiliauskas, T. lakimov, A. Zakharov, 
M. Syvajarvi, and R. Yakimova, Carbon 57, 477 (2013). 
Y. Lin, A. Valdes- Garcia, S. Han, D. Farmer, L Meric, 
Y. Sun, Y. Wu, C. Dimitrakopoulos, A. Grill, P. Avouris, 
et al.. Science 332, 1294 (2011). 

S. Hertel, D. Waldmann, J. Jobst, A. Albert, M. Al- 
brecht, S. Reshanov, A. Schoner, M. Krieger, and H. We- 
ber, Nature Comm. 3, 957 (2012). 

C. Berger, Z. Song, T. Li, X. Li, A. Ogbazghi, R. Feng, 
Z. Dai, A. Marchenkov, E. Conrad, N. Phillip, et al., J. 
Phys. Chem. B 108, 19912 (2004). 

P. First, W. De Heer, T. Seyller, C. Berger, J. Stroscio, 
and J. Moon, MRS Bulletin 35, 296 (2010). 
S.-H. Ji, J. B. Hannon, R. M. Tromp, V. Perebeinos, 
J. Tersoff, and F. M. Ross, Nature Materials (2011). 
K. V. Emtsev, F. Speck, T. Seyller, L. Ley, and J. D. 
Riley, Phys. Rev. B 77, 155303 (2008). 
F. Ming and A. Zangwill, Phys. Rev. B 84, 115459 
(2011). 

K. Renter and M. Scheffler, Phys. Rev. B 65, 035406 
(2001). 

R. M. Tromp and J. B. Hannon, Phys. Rev. Lett. 102, 
106104 (2009). 
C. Riedl, U. Starke, J. Bernhardt, M. Franke, and 



[is: 

[19 

[20; 

[21 

[22; 
[23; 

[24 
[25; 

[26; 

[27 



[28; 
[29; 

[30 
[31 



[32' 
[33^ 

[34' 
[35' 
[36' 



[37 
[38; 
[39; 
[40; 
[41 
[42; 

[43; 
[44; 



K. Heinz, Phys. Rev. B 76, 245406 (2007). 
J. B. Hannon (2012), private communication. 

A. Tkatchenko and M. Scheffler, Phys. Rev. Lett. 102, 
073005 (2009). 

J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 

Lett. 78, 1396 (1997). 

S. Grimme, Comput. Mol. Sci. 1, 211 (2011). 

M. Dion, H. Rydberg, E. Schroder, D. C. Langreth, and 

B. L Lundqvist, Phys. Rev. Lett. 92, 246401 (2004). 

K. Lee, E. D. Murray, L. Kong, B. L Lundqvist, and 
D. C. Langreth, Phys. Rev. B 82, 081101 (2010). 
J. Klimes, D. R. Bowler, and A. Michaelides, J. Phys.: 
Condens. Matter 22, 022201 (2010). 
G. X. Zhang, A. Tkatchenko, J. Paier, H. Appel, and 
M. Scheffler, Phys. Rev. Lett. 107, 245501 (2011). 
A. Tkatchenko, R. A. DiStasio, R. Car, and M. Scheffler, 
Phys. Rev. Lett. 108, 236402 (2012). 
See supplementary material. The supplementary material 
includes reference data for the numerical convergence of 
our calculations, calculated lattice parameters and other 
cohesive properties of the bulk reference phases 3C-SiC, 
diamond Si, diamond and graphite C, without and with 
zero-point corrections for the functionals used here, and 
structural details of the exemplary defects as introduced 
in Ref. ;58]. 

V. Blum, R. Gehrke, F. Hanke, P. Havu, V. Havu, 
X. Ren, K. Renter, and M. Scheffler, Comp. Phys. Com- 
mun. 180, 2175 (2009). 

V. Havu, V. Blum, P. Havu, and M. Scheffler, J. Comp. 
Phys. 228, 8367 (2009). 

T. Auckenthaler, V. Blum, H. Bungartz, T. Huckle, 
R. Johanni, L. Kramer, B. Lang, H. Lederer, and 
P. Willems, Parallel Computing 37, 783 (2011). 
U. Starke, J. Schardt, J. Bernhardt, M. Franke, 
K. Renter, H. Wedler, K. Heinz, J. Furthmiiller, 
P. Kackell, and F. Bechstedt, Phys. Rev. Lett. 80, 758 
(1998). 

Surf. Sci. 215, 111 (1989). 

U. Starke, J. Schardt, J. Bernhardt, M. Franke, and 
K. Heinz, Phys. Rev. Lett. 82, 2107 (1999). 
Z. Li and R. C. Bradt, J. Mat. Sci. 21, 4366 (1986). 
Y. Baskin and L. Meyer, Phys. Rev. 100, 544 (1955). 
T. Seyller, A. Bostwick, K. Emtsev, K. Horn, L. Ley, 
J. McChesney, T. Ohta, J. Riley, E. Rotenberg, and 
F. Speck, phys. stat. sol. (b) 245, 1436 (2008). 
J. B. Hannon, M. Copel, and R. M. Tromp, Phys. Rev. 
Lett. 107, 166101 (2011). 

T. Ohta, A. Bostwick, T. Seyller, K. Horn, and E. Roten- 
berg, Science 313, 951 (2006). 

S. Kim, J. Ihm, H. J. Choi, and Y. W. Son, Phys. Rev. 
Lett. 100, 176802 (2008). 

F. Varchon, P. Mallet, J.-Y. Veuillen, and L. Magaud, 
Phys. Rev. B 77, 235412 (2008). 

G. Sclauzero and A. Pasquarello, Phys. Rev. B 85, 
161405(R) (2012). 

W. Chen, H. Xu, L. Liu, X. Gao, D. Qi, G. Peng, S. C. 
Tan, Y. Feng, K. P. Loh, and A. T. S. Wee, Surf. Sci. 
596, 176 (2005). 

W. De Heer, T. Seyller, C. Berger, J. Stroscio, and 
J. Moon, MRS Bulletin 35, 296 (2010). 

C. Berger, J. Veuillen, L. Magaud, P. Mallet, V. Ole- 
vano, M. Orlita, P. Plochocka, C. Faugeras, G. Martinez, 
M. Potemski, et al.. Int. J. Nanotechnol. 7, 383 (2010). 



[45] S. Goler, C. Coletti, V. Piazza, P. Pingue, F. Colangelo, 

V. Pellegrini, K. V. Emtsev, S. Forti, U. Starke, F. Bel- 
tram, et al.. Carbon 51, 249 (2013). 
[46] L. H. de Lima, A. de Siervo, R. Landers, G. A. Viana, 

A. M. B. Goncalves, R. G. Lacerda, and P. Haberle, Phys. 

Rev. B 87, 081403(R) (2013). 
[47] W. Norimatsu and M. Kusunoki, Chem. Phys. Lett. 468, 

52 (2009). 
[48] F. Varchon, R. Feng, J. Hass, X. Li, B. Nguyen, C. Naud, 

P. Mallet, J. Veuillen, C. Berger, E. Conrad, et al., Phys. 

Rev. Lett. 99, 126805 (2007). 
[49] W. Lu, J. J. Boeckl, and W. C. Mitchell, J. Phys. D: 

Appl. Phys. 43, 374004 (2010). 
[50] J. Roehrl, M. Hundhausen, K. Emtsev, T. Seyller, 

R. Graupner, and L. Ley, Appl. Phys. Lett. 92, 201918 

(2008). 
[51] R. Berman and F. Simon, Z. f. Elektrochem.: Ber. Bun- 

senges. Phys. Chem. 59, 333 (1955). 
[52] M. T. Yin and M. L. Cohen, Phys. Rev. B 29, 6996 



(1984). 
[53] F. Lazarevic, L. Nemec, V. Blum, and M. Scheffler, to be 

published (2013). 
[54] G. Sclauzero and A. Pasquarello, Diamond Rel. Materials 

23, 178 (2012). 
[55] A. Mattausch and O. Pankratov, Phys. Rev. Lett. 99, 

076802 (2007). 
[56] O. Pankratov, S. Hensel, and M. Bockstedte, Phys. Rev. 

B 82, 121416 (2010). 
[57] U. Starke and C. Riedl, J. Phys.: Condens. Matter 21, 

134016 (2009). 
[58] Y. Qi, S. H. Rhim, G. F. Sun, M. Weinert, and L. Li, 

Phys. Rev. Lett. 105, 085502 (2010). 
[59] P. Sutter, Nature Materials 8, 171 (2009). 
[60] J. Hass, R. Feng, J. E. Millan-Otoya, X. Li, M. Sprinkle, 

P. N. First, W. A. de Heer, E. H. Conrad, and C. Berger, 

Phys. Rev. B 75, 214109 (2007). 



